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We report on the occurrence of multiple hopping and retrapping of a Brownian particle in a tilted 
washboard potential. The escape dynamic has been studied experimentaUy by measuring the switch- 
ing current distributions as a function of temperature in a moderately damped NbN/MgO/NbN 
Josephson junction. At low temperatures the second moment of the distribution increases in agree- 
ment with calculations based on Kramers thermal activation regime. After a turn-over temperature 
T* , the shape of the distributions starts changing and width decreases with temperature. We an- 
alyze the data through fit of the switching probability and Monte Carlo simulations and we find a 
good agreement with a model based on a multiple retrapping process. 

PACS numbers: 05.40.Jc, 74.50.-|-r, 85.25.Cp 



INTRODUCTION 



Research on superconducting quantum sys- 
tems with potentials for qubits applications has 
boosted interest on several complementary as- 
pects of coherence and dissipation. There is a 
growing evidence of the occurrence of a mod- 
erately damped regime (MDR) in supercon- 
ducting Josephson junctions (JJs) of various 
materials^"— . If we use the Q = uipRC parame- 
ter as a measure of dissipation in a junctionii^, 
where Up = (2e/c/?iC)^/^ represents the plasma 
frequency at zero bias current, Ic is the junction 
critical current and C and R are the junction 
capacitance and resistance, a MDR is present 
for 1 < Q < 5. This regime is quite distinct 
from the well-known case of underdamped sys- 
tems (Q > 10)^'^", and apparently quite com- 
mon in junctions characterized by low Ic- In 
view of a more and more relevant use of nan- 
otechnologies in quantum superconducting elec- 



tronics and therefore of low values of Ic, studies 
on MDR can offer novel insights on dissipative 
effects on Josephson junctions, and inspire ap- 
propriate designs to respond to specific circuit 
requirements. 

Our analysis is based on measurements on low 
critical current density (Jc) NbN/MgO/NbN 
junctions. NbN is a material of great interest for 
sensor applications, as documented by several 
works both on junctions and nanowiresii"— , 
and it guarantees both fast non equilibrium 
electron-phonon relaxation times r < lOps and 
higher gap values, when compared with tradi- 
tional junction technologies based on Nb, Al 
and Pbi^. Low-Jc NbN/MgO/NbN devices 
may contribute to set a more comprehensive 
NbN platform, and constitute a non-trivial ex- 
tension to thicker barriers of the more estab- 
lished high- Jc NbN junctions, usually designed 
for superconducting digital circuit a^^d"^ . 

In this work we demonstrate the occurrence 
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of a phase diffusion regime (PDR) induced by 
low- Jc in large area junctions. Areas are one 
or two order of magnitudes larger than those 
of junctions where PDR has been previously 
observed^'^, and Josephson (Ej) and Coulomb 
(Ec) energies entering in the tunneling process 
will scale with size in a different manner because 
of the diverse Jc and specific capacitance (Cs) 
values. PDR prevails over thermal activation 
at temperatures above a threshold T* of about 
1.5K . Experimental results are very close to 
theoretical predictionsi*^ such to provide a re- 
liable estimation of Q = 2.7. We compare our 
results with the model of Fenton and Warburton 
(FW)^ which condenses ideas on phase diffusion 
of the last 20 years, offering a reliable method- 
ology to evaluate levels of dissipation in MDR, 
in analogy with what well established for under- 
damped junctions^ii^. In particular FW model 
provides additional criteria to study low Q junc- 
tions to prove their MDR nature based on the 
asymmetry of the switching current histograms 
or its skewness 7, i.e. the ratio m^/a^ where 
777.3 is the third central moment of the distribu- 
tion, which are closely followed by experimen- 
tal data. A physical picture emerges of mod- 
erately damped junctions, with a damping sub- 
stantially independent of the frequency and able 
to sustain macroscopic quantum tunnelingiS. at 
lower temperatures. 



II. THEORETICAL BACKGROUND 

The significance of the analysis of phase dif- 
fusion phenomena extends to the more general 
problem of the motion of a Brownian parti- 
cle in a periodic potential, which is behind the 
Resistively Shunted Junction (RSJ) model for 
the dynamics of the phase difference across a 
Josephson junction^'^. It can describe many 
different physical phenomena, i.e. transport on 
crystalline surfaces-- , rotating dipoles in exter- 
nal fielda^S, charge density wavea^, and particle 
separation by electrophoresis^S. Kramers^'^ first 
studied the one-dimensional escape dynamic of 
a quantum Brownian particle in presence of 
damping, in the low and high damping limits. 




FIG. 1: Tilted periodic potential and the equivalent 
circuit of the RSJ model. 

and subjected to a tilted potential U{z) 

U{z)^V{z)-Fz (1) 

where V{z) — V{z + P) denotes a periodic po- 
tential of period P and F is an external static 
force. The dynamics of the phase difference <p 
of the superconducting order parameter across 
the junction of a current-biased JJ is equivalent 
to the Brownian motion of a particle of " mass" 
C (the junction capacitance) in a tilted periodic 
potentialii^ 

U{^) = -Ej(^os^+j^ (2) 

In this case the external static force is con- 
trolled by the bias current I and is given by 
Ej{I/Ic) where Ej — hic/'i-e is the Josephson 
energy. The motion of the particle along the po- 
tential is also subject to friction, whose strength 
can be characterized by the above mentioned 
junction quality factor Q. In this model the su- 
perconducting branch of the junction current- 
voltage (I-V) characteristic corresponds to the 
confinement of the particle in one well of the po- 
tential. The escape from this metastable state 
corresponds to the appearance of a finite volt- 
age across the junction. As it is shown in fig. 
[U in case of low damping the escaped particle 
gains sufficient energy to roll down the potential 
in the so-called running state, meanwhile if the 
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damping is sufficiently high, escape due to ther- 
mal hopping does not necessarily lead to runway 
down the tilted potential^i^. Following an event 
of escape the particle may travel down the po- 
tential for a few wells and then be retrapped in 
one of the following minima of the potential^. 
At low bias the process of escape and retrapping 
may occur multiple times generating extensive 
diffusion of the phase until an increase of the 
tilt of the potential, due to a change in the bias 
current, raises the velocity of the particle and 
the junction can switch to the running statc'\ 

This phenomenon of phase diffusion is more 
evident when studying the temperature depen- 
dance of the switching probability^^: 



T(I) 



nn 
di'/dt 



di' 



(3) 



where the rate of escape due to thermal activa- 
tion is given by Kramers formulai^^ 



TA 



(I) 



at— — exp 



^""') (4) 



being AJ7 the height of the energy barrier be- 
tween consecutive potential wells and the pref- 
actor is at = 4/[(l + QhT/l.SAU)^^^ + l]^. 
For underdaniped junctions the width a of the 
switching distributions has a monotone depen- 
dence on temperature as cr oc T^^^. In the 
case of moderately damped junctions, due to 
the retrapping process the switching dynamic 
is modified; below a turn-over temperature T* 
the width of the switching distributions follows 
the usual a oc T^/^ , while for T > T*, cr is 
reduced with increasing temperature. In this 
paper we also report on the observation of this 
behavior which we fit using the FW model. 



III. SAMPLES FABRICATION 
PROCESS AND EXPERIMENTAL 
SETUP 

In NbN/MgO/NbN JJs a 1 nm thick barrier 
provides Jc of about SA/cm^, which is the low- 
est value ever reported for NbN based junctions. 
For circular junctions with a diameter of 10/im, 



Ic is about 1-2 and falls under the crite- 
ria of the moderately damped regime, as exten- 
sively discussed below. In the trilayer, epitaxi- 
ally grown at ambient temperature on a single- 
crystal MgO substrate^S, the NbN base (BE) 
and counter electrode (CE) are both 200nm 
thick and were deposited using DC magnetron 
sputtering with a Nb target in a mixture of 5 
parts argon and 1 part nitrogen gas. The MgO 
barrier is about l.Onm thick and was deposited 
by rf sputtering. This step was followed by a 
reactive ion etch (RIE) for junction definition 
and by the deposition of a MgO insulating layer 
patterned by a lift-off process. The process is 
concluded with the deposition of a 350nm NbN 
wiring layer which was patterned and defined 
by RIE. The realized junctions have a supercon- 
ducting transition temperature of about 16.6 K 
for both electrodes. A more detailed descrip- 
tion of the fabrication process can be found in 
a paper from Kawakami et al.— . 

The junctions have been tested through mea- 
surements of the I-V characteristics (see fig. 
[2])and of switching current distribution (see next 
section). From the magnetic field dependence 
of the critical current (shown in the inset of 
fig. m we estimated the London penetration 
depth at 300mK to be about Xl = 190nm which 
is in good agreement with previously measured 
values for epitaxially grown NbN27'2^ therefore 
for a barrier of thickness t=l nm the Joseph- 
son penetration depth turns out to be^ Xj = 
{hcySneJcdy^^ = 150/im, where d = 2Xl + 1. 

To study the escape rates of a 
NbN/MgO/NbN Josephson junction we 
have thermally anchored the sample to the 
mixing chamber of a He3/He4 Oxford dilution 
refrigerator and performed measurements of the 
junction switching current probability. Figure 
[3] shows a block diagram of the experimental 
setup including room temperature electronic 
and filtering. The room temperature circuits 
have been optimized to minimize the effect 
of unwanted noise. In order to avoid ground 
loops and noise pick up, the whole experiment 
is designed to have a single groun d^^'^" . All 
connections to the chip are floating and are 
only capacitively coupled to ground through 
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FIG. 2: Current- Voltage characteristic of a Wfim 
diameter JJ measured at 290mK. The inset shows 
the magnetic field dependence of the Josephson tun- 
nel current, with magnetic field values normalized 
to the first minimum. 



the filters. All grounded signal sources pass 
through battery powered unity gain isolation 
amplifiers that effectively disconnects this 
signal from the earth ground. The current 
paths are all designed to be symmetric with 
respect to the chip to reduce the effect of com- 
mon mode noisc"^^. This allows for sufficient 
decoupling from the ground while keeping the 
amplifier from saturating due to charging by 
providing a return current path. The amplifiers 
are designed to have 100 Mil of resistance 
between their inputs and the common of the 
battery circuit. All signals entering the fridge 
are isolated, shielded and filtered allowing the 
dilution refrigerator itself to act as an rf-shield 
for the cold portion of the experiment. In our 
system we use a room temperature electro- 
magnetic interference filter stage followed by 
low pass RC filters with a cut-off frequency of 
1.6MHz anchored at 1.5 K. Further filtering is 
provided by a combination of copper powder— 
and twisted pair filters^ thermally anchored at 
the mixing chamber of the dilution refrigerator. 



IV. DATA AND ANALYSIS 

The signal sequence used to measure the 
switching current distribution (SCD) is shown 
in the inset of fig HI The bias current of the 
junction is ramped at a constant sweep rate 
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FIG. 3: (color online) Schematic of measurement 
electronics including thermal anchoring and the var- 
ious stages of filtering. 



dl/dt = 122nA/s, the voltage is measured using 
a low noise differential amplifier and is fed into 
a threshold detector which is set to generate a 
pulse signal when the junction switches from 
the superconducting state to the finite volt- 
age state. This signal is used to trigger a fast 
volt meter to record the value of the switching 
current^'*. This procedure is repeated at least 
10'' times at each temperature, which allows to 
compile a histogram of the switching currents. 
In figlH we report the SCD curves collected 
over a wide range of temperatures in absence 
of an externally applied magnetic field. Dis- 
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FIG. 4: Switching current probability distribution 
at B = OG for different bath temperatures. The 
inset shows the signal sequence used to acquire the 
data. 



tinctive fingerprints of phase diffusion, due to 
multiple hopping and retrapping, can be found 
in the temperature dependence of the width a 
of the SCD curves, which is shown in Fig. [5l 
The most striking effect observable in FiglS] is 
the appearance of an anticorrelation between 
the temperature and the width of the switch- 
ing distributions^r— . At low temperatures the 
a follow the expected T'^/'^ dependence, devi- 
ations are evident in proximity and above a 
"critical temperature" T* where the tempera- 
ture derivative of cr(T) becomes negative. Ex- 
perimental data (upper frame in FiglSJ are well 
reproduced by the expected values (lower frame 
in FiglSJ , calculated on the basis of the physical 
arguments of phase diffusion. 

The simulations shown in Fig (5] in particular, 
are based on the recent work on phase diffusion 
by Fenton and Warburton^. The phase differ- 
ence ip{t) is a solution of the following Langevin 
differential equation: 



(5) 



Times t are normalized to l/w^; u is the bias 
current normalized to critical current Ico and t n 
is a Gaussian correlated thermal noise current, 
i.e.: 



{iN{t), iN{t')) = ^2TTkBT/Qho^o5[t - t'). 



Stochastic dynamics is simulated by integrat- 
ing the above Langevin equation by a Bulirsh- 
Stoer integrator using as noise generator the 
cernlib routine RANLUX^. Simulations have 
been carried out for different temperatures and 
dissipation values. 

The multiplicity of switching modes between 
the running and the trapped states raises a 
problem of how to define an escape event. In our 
simulations the condition to define the switch is 
yii-iT) > V{l,0)/2, where V represents the av- 
erage velocity of the phase-particle in the wash- 
board potential. In other words the particle 
spends in the running state more than 50% of 
observation time. Typical runs for simulations 
of Eq.® will last from 4 • 10*^ to 6 • 10^ nor- 
malized time units, i.e., 6 • 10^ to 9 • 10^ plasma 
periods. Observation time for each point gen- 
erated in the TV characteristics is 2 • 10'* time 
units, which is a long enough time to ensure that 
the average time spent in running/zero voltage 
state does not vary as a function of the ob- 
servation time*. To obtain the SCD we have 
simulated a number of escape events between 
3000 and 5000, which is similar to the number 
of counts experimentally measured. Simulated 
curves of cr vs T for different values of the mag- 
netic field are plotted in Fig. 5. The magnetic 
field works as a knob to tune T* and provides 
an additional validity test for the estimate of 
Q = 2.7 ± 0.1. From fitting of the SCD and 
of the moments of the distributions at temper- 
atures below T* we have estimated the value 
of the zero temperature critical current to be 
Ico — 1-91 ± 0.03/iA. We have also estimated 
values for the junction capacitance and plasma 
frequency and obtained C = O.SpF and ujp ~ 
22 GHai^. 

The quality of the fitting procedure is even 
more significant if we consider that we do not 
have any degree of freedom associated to a pos- 
sible frequency dependence of Q as occurring in 
other experiments^'^'^. For instance the proce- 
dure used for our analysis is different from that 
used by Mannik et ali^, where retrapping prob- 
ability is calculated independently for different 
dissipations. An analysis of full escape rate is 
made by combining the V escape probabil- 
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FIG. 5: Top frame: temperature dependence of the 
standard deviation, ct, of tiie switching distributions 
for B = OG (squares), B = 3.05G (circles), and 
B — 6.09G (triangles). Bottom frame: a numerical 
simulation of the data. Vertical dotted lines have 
been inserted in correspondence of the values of T* 
including error bars. Data and numerical simula- 
tions are in good agreement in the whole tempera- 
ture range and for all magnetic fields within error 
bars.. 




0.1 1 
Temperature (K) 

FIG. 6: Temperature dependence of mean cur- 
rent, /, of the switching distributions for B = OG 
(squares), B = 3.05G (circles), and B = 6.09G (tri- 
angles). The solid lines are the calculated values 
of I considering thermal activation and no retrap- 
ping. The inset shows the value of T* vs the applied 
magnetic field. 



ity with retrapping probability which allows to 
extract the resistance by fitting the escape rate 
curves. 

In Fig. [6] the mean values / of the SCD (data 
points) are plotted along with the expected val- 
ues (solid lines) without taking into account re- 
trapping effects. Due to the onset of retrap- 
ping events it is necessary to provide a larger 
tilt to the energy potential to allow the system 
to switch to the running state. Discrepancies 
at higher T demonstrate that the experimen- 
tal values of / at higher temperatures (above 
T*) are greater than the predicted values, which 
only consider the effects of thermal activation. 
Due to the dependence of T* with the external 
magnetic field, as shown in the inset of Fig. |6l 
this latter effect is particularly evident for the 
data taken in presence of a 6.09G static mag- 
netic field, where the onset of retrapping occurs 
at a lower temperature. 

Phase diffusion also appears in the escape 
rates F, shown in figure[7]as a function of the ra- 
tio between the barrier height and the thermal 
energy. The escape rates are calculated from 
the switching distributions using eq. [3] In the 
thermal activation regime the distributions are 
asymmetric and skewed to the left, and F val- 
ues all fall onto the same line, as it is the case 
for the reported data from T=0.3K to 1.56K. 
Retrapping processes cause a progressive sym- 
metrization of the switching distribution, as it 
can be seen from the inset in the bottom left 
corner of Fig. [71 and a bending in the F vs 
u = AU/kBT. We use the same procedure pre- 
viously described to evaluate the numerically 
simulated escape rates F as function of reduced 
barrier height u. Numerical data have been ob- 
tained by a polynomial fit of numerical escape 
rates in order to compare it with experiments. 
The same value of the Q factor is obtained by 
fitting the T{u) curves, shown in dashed lines in 

Fig.m 

The symmetrization of the switching distri- 
bution due to the interplay between escape and 
retrapping events can be clearly observed by 
plotting, as a function of temperature, the skew- 
ness of the distributions 7^^. For the data with 
no external magnetic field, we report such plot 
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FIG. 7; Escape rates (symbols) as a function of the 
barrier height at zero magnetic field, for tempera- 
tures near T* , and numerical simulations (dashed 
lines). Below T* the fit has been calculated using 
Kramers formula for thermal activation (solid line). 
The inset shows the experimental value of the skew- 
ness of the switching distribution. 

in the inset of fig. [71 For the lowest tem- 
peratures we obtain 7 = — 1, which is consis- 
tent with the case of switching current distri- 
butions in the quantum or thermal regime. As 
the temperature increases the distributions be- 
come more and more symmetric as 7 tends to 
zero. It should be noted that for these data the 
temperature T* at which the derivative of 'j[T) 
changes sign is equal to 1.62 ± 0.03 K and that 
the skewness starts increasing already at about 
\.2K, which is a clear indication that the on- 
set of retrapping phenomena occurs well below 



V. DISCUSSION AND CONCLUDING 
REMARKS 

A change in the sign of the derivative of 
the second moment of the distribution at a 
turn-over temperature T* and a modification 
of the shape of the distributions at tempera- 
ture around T* are robust signatures of the 
phase diffusion regime, and also occur in our 
NbN/MgO/NbN junctions, as discussed in the 
previous section. A non exhaustive list of 
Josephson devices (including ours) that have 



displayed a similar PDR behavior is reported 
in table U along with the most relevant device 
parameters. 

The papers on which tableUis based report on 
similar experimental results but their interpre- 
tation differs in few assumptions, as properly 
pointed by Fenton and Warburtor*^ . For in- 
stance Kivioja et al.^ interpreted their results 
within the semiclassical model of Larkin and 
Ovchinniko v^^i'^^ Since in dc-squids there are few 
energy levels and the hypothesis of continu- 
ous energy spectrum is not valid, they used a 
model which takes account both phase diffu- 
sion and level quantization. On the other hand 
this model, which assumes separated levels in 
the metastable well, is not properly valid for 
a single Josephson junction since the number 
of energy levels is large and the separation is 
smaller than their width. Mannik et al.'^ and 
Bae et al<^ calculated the retrapping probabil- 
ity through Monte Carlo simulations and in- 
cluded frequency dependent damping. The au- 
thors expressed the net escape rate as a sum 
of probabilities of multiple escape-retrapping 
events based on thermal escape rate and re- 
trapping probability. The probability of re- 
trapping is considered as a time-independent 
quantity which is in contrast with the work of 
Ben- Jacob et al<^ in which retrapping is mod- 
eled by a rate and therefore by a probability 
increasing proportionally to the time spent in 
the running state. The very good fitting of ex- 
perimental curves obtained in this work using 
the FW approacli^ confirms the occurrence of a 
multiple-retrapping regime with a large number 
of escapes of duration of F^^ , and in particular 
it confirms that the time dependence of the re- 
trapping probability cannot be ignored. If the 
fast scattering time r plays a role in diffusive 
process is a topic of further investigations. 

As it can be seen from table HI our experiment 
confirms that, independently of the physical size 
of the device, all the junctions exhibiting phase 
diffusion over a large range of materials and ge- 
ometry have a low critical current, 2 < Q < 5 
and 12 < AU(l)/kBT* < 18, which are there- 
fore the relevant parameters signaling the in- 
surgence of multiple escape and retrapping in a 
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TABLE I: Comparison of device parameters. 



Author 


Device Structure 


R [Vt) 


Ico {^J^.A) 


Area {^.m^) Q{I = 0) 


AU{I)/kBT* 


Ej/Ec 


Tliis work 


NbN/MgO/NbN J J 


65 


1.91 


78.5 


2.7 


17 


14760 


Kivioja et al? 


Al/AlOx/Al dc SQUID 


500 


0.2 


1.0 


3.9 


14 


515 


Kivioja et al^ 


Al/AlOx/Al JJ 


230 


0.63 


2.6 


3.6 


18 


2110 


Mdnnik et al^ 


Nb/AlOx/Nb dc SQUID 


70 


4.25 


1.0 


2.4 


15 


9850 


Mdnnik et al^ 


Nb/AlOx/Nb dc SQUID 


70 


2.9 


2.9 


3.3 


17 


19420 


Bae et al^ " 


Bi-2212 Intrinsic JJ 


62 


1.26 


7.3 


2.2 


14 


10710 


Yu et al^ 


Nb/AlOx/Nb JJ 


1800 


0.122 


0.4 


4.8 


NA 


62 


Yu et al}° 


Nb/AlOx/Nb JJ 


315 


0.48 


1.5 


3.3 


12 


950 



"In tills paper the authors estimated the fit parameters 
to be temperature dependent. Here we report the values 
at the lowest experimental temperature T=1.5K. 



washboard potential. 

The possibihty to have extremely low critical 
current density can be functional to investigate 
phase dynamics at extreme conditions. An ex- 
ample is given by a recent experiment on sub- 
micron Nb/AlOx/Nb junctions^. Data show 
an anomalous (t(T) dependence with a negative 
da/dT over the entire temperature range and 
a saturation at low temperatures. This regime 
can be achieved by engineering junctions with 
low critical current, such that the turn over tem- 
perature T* is lower or comparable to the quan- 
tum crossover temperature Tcrosi^- In this case 
the enhancement of a when increasing the tem- 
perature, characteristic of the thermal regime, 
is not observed. Junctions with intrinsically low 
critical current density, such as the one reported 
in the present work, could represent an interest- 
ing term of comparison to study these kind of 
unconventional regimes using standard microm- 
eter junctions. 

In conclusions we have proved that low 
NbN/MgO/NbN JJs are characterized by a 
transition from thermal activation regime to 
phase diffusion. This is consistent with what 



has been observed in other types of junctions 
with similar values of Q and A[/(/)/fc_BT*. The 
experimental results are well described by a nu- 
merical model involving a frequency indepen- 
dent damping which demonstrates an efficient 
way to estimate the dissipation in moderately 
damped JJs. 
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